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ABSTRACT 

Context. The surface brightness fluctuation (SBF) method is a powerful tool to derive distances to galaxies for which single stars cannot be 
resolved. Up to now, the method has been calibrated mainly at red colours due to the intrinsic faintness of blue early-type galaxies. 
Aims. In this paper we address the /-band calibration of the SBF method at blue colours, the regime of dwarf elliptical galaxies (dEs). 
Methods. We present deep and wide-field VI photometry of the central Fornax cluster obtained at Las Campanas Observatory. With these data 
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we perform an SBF analysis of 25 dEs in the range -16.5 < My < -1 1.2 mag, 0.8 < {V - I)q < 1.10 mag. Our colour calibration is accurate to 
better than 2%. For the calibration analysis we exclude eight dEs whose SBF measurement was affected by poor seeing (FWHM =^ 1"). 
Results. Our SBF data are inconsistent at the 3cr level with a colour independent absolute SBF magnitude M/, presenting a problem 
• for SBF models that predict such a "flat" relation. There is a weak indication in our data (1.8cr) to favour a two-branch calibration over 
a one-branch calibration with broad scatter in the Mj - (V - I) plane. We obtain the following one-branch empirical SBF calibration: 
Q , M/ = -2.13(±0.17) + 2.44(±1.94) x [(V - /)o - 1.00] mag. We deduce a 0.34 ± 0.14 mag cosmic scatter of M/, which is significantly larger 
5_j • than found at redder colours. This is in agreement with those theoretical SBF models that predict Mj to be more sensitive to age-metallicity 
variations in the blue than in the red. We find evidence that the fainter galaxies in our samples contain younger and more metal-rich stellar 
d populations than the brighter ones. The application of our empirical calibration to published SBF measurements of Hydra and Centaurus 
ILJ cluster dEs leaves the distances to both clusters unchanged. 



^ Key words, galaxies: clusters: individual: Fornax cluster - galaxies: dwarf - galaxies: distances and redshift - techniques: photometric 



1. Introduction 

1.1. The principle of SBF 

The measurement of surface brightness fluctuations (SBF) in a 
galaxy image was first suggested by Tonry & Schneider ( 1988 1 
as a useful tool to measure distances to galaxies which are not 
resolved into single stars (see also Jacoby et al. il992i for a re- 
view). The SBF method is based on two facts: 

A) the number of stars contained in each resolution element of 
a galaxy image is finite. 

B) the number of stars contained in each resolution element is 
on average constant along one isophote. 

From these two facts, a third fact is deduced: the number of 
stars per resolution element along one isophote has a statistical 
fluctuation equal to the square root of the average star num- 
ber. The amplitude of those surface brightness fluctuations nor- 
malised to the underlying mean galaxy luminosity is inversely 
proportional to distance and can therefore be used as a distance 
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indicator. This normalised amplitude converted to magnitude 
is in the framework of SBF investigations referred to as appar- 
ent fluctuation magnitude m. It corresponds to the luminosity 
weighted mean luminosity of stars in the observed stellar pop- 
ulation (see e.g. Tonry & Schneider 1988 and Liu et al. l2000t . 
The highest luminosity stars in an old stellar population are the 
red giants. Consequently, the average number of stars per unit 
surface brightness is lowest in the colour regime of red giants. 
SBF measurements have therefore mainly been performed in 
red optical pass-bands (e.g. Tonry & Schneider [198 8 Bothun 
et al. 1991 Tonry et al. flWTI Jerjen et al. 2001 Meiet al.|2O05} 
and in the near infrared (e.g. Luppino & Tonrv [T993l Pahre & 
Mould EHHI Liu & Graham 2001 Liu et al. IMHl Jensen et 
al. l20n^ . bofli from the ground and with HST. 

1.2. SBF calibrations 

Of course, any distance determination method needs to be cali- 
brated in order to derive absolute distances. In their pioneering 
work, Tonry et al. 119971 1200 U have carried out an extensive 
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survey to measure the apparent /-band SBF magnitude m/ 
for several hundred bright elUpticals and bulges of spiral 
galaxies in 22 nearby galaxy groups and clusters within 40 
Mpc. Using distances derived from Cepheids or SN la in the 
respective galaxies or the clusters they are associated with, 
they obtained a purely empirical relation between the absolute 
I-band fluctuation magnitude M/ and the dereddened colour 
(V - /)o, determined in the range 1.05 < (V - /)o < 1.3: 

M, = -1.74 + 4.5 x((y-/)o- 1.15) mag (1) 

The sign of this slope is in agreement with the expectation that 
stars in blue stellar populations have higher average luminosi- 
ties than stars in red populations. The value of the slope agrees 
reasonably well with theoretical expectations, while zero-point 
predictions from SBF models scatter significantly around the 
value Tonry obtained (see also Sect. |5] and the discussion in 
Tonry et al. 2001 1. A limitation of the Tonry calibration is that 
for observational reasons it is restricted to rather luminous and 
large galaxies. It does not include any dwarf galaxies and does 
therefore exclude the blue colour range. The advantage of a 
calibration for bluer colours is the high number of dwarf galax- 
ies as compared to the sparser red giant galaxies (e.g. Karick 
et al.|20p3 Peng et al. 12^051 Hilker et al."2003^, increasing a 
lot the number of accessible independent distance estimators. 
This gain of measurement points is partially compensated by 
the faint surface brightness and smaller intrinsic size of dEs, 
which lowers the measurement accuracy (see Sect. 13.2.11 . 
In Mieske & Hilker (2003b i and Mieske et al. (2005 1 we de- 
rived /-band SBF distances to the Centaurus and Hydra clusters 
and discussed how the Great Attractor hypothesis (e.g. Lynden- 
Bell et al. fl988l Tonry et al. EOOOl fits into the derived pecu- 
liar velocities. These measurements were mainly based on dis- 
tances to dEs with colours (V - /)o < 1.10 mag, i.e. close to 
or bluer than the blue limit of the Tonry calibration. A 50% 
shallower calibration relation than the Tonry one was assumed 
for the bluest dEs, inspired by the scatter of theoretical models 
in that colour range (see also Sect. The main result of our 
studies was that the Hydra cluster is somewhat closer to us or 
at most at the same distance as Centaurus. The coiTespondingly 
low Centaurus and high Hydra peculiar velocity was consistent 
with a Great Attractor in close projection to and slightly be- 
hind Hydra. An empirical or semi-empirical SBF calibration in 
the blue regime would be very helpful to improve the accuracy 
of the derived distances, especially the relative one between 
Hydra and Centaurus. 

Jerjen et al. ( 1998 2000 2(50?) derived a semi-empirical SBF 
calibration between Mr and (B - R) in the range 1.0 < {B - 
R) < 1.3 mag by targeting nearby dwarf elliptical galaxies 
with known distances. These authors based their calibration on 
the theoretical predictions of Worthey ( 1994 1 models and used 
their observational data to adjust the zero-point of those mod- 
els. In Jerjen et al. (2001 1 they found that the models predicted 
too faint SBF amplitudes by 0.13 mag. Furthermore, they pro- 
pose two separate branches in the Mr-{B - R\) plane at blue 
colours. At a given colour, becomes brighter with younger 
and more metal-rich stellar population. This means that the 
age-metallicity degeneracy of Mr is partially broken at blue 



colours. 

In the course of the ACS Virgo Cluster survey, Mei et al. (I2005> 
found a significant flattening of the SBF-colour relation in M, 
vs. (g - z) for \.Q < ig - z) < 1.3. They obtained a significant 
slope in that colour range (4.5cr), but did not find evidence for 
a bimodality of M^ values. 

In the infrared, Jensen et al. (I2003> presented an HST based 
SBF calibration between MF160W and (V - /)o in the range 
1.05 < (y - /)o < 1.24 mag, which was based only on giant 
galaxies and yielded a comparable slope to that in equation Q. 
For ground based infrared bands, Jensen et al. J1998> and Liu 
et al. ( 2002 1 presented calibrations between K and {V -!)(), also 
based purely on giant galaxies. The small sample of Jensen et 
al. in the range 1.15 < (y-/)o < 1.27 mag was consistent with 
a zero slope between (V - /)(> and Mk- Liu et al. obtained a 
significant slope of 3.6 + 0.8 between both entities in the range 
1.05 <(y-/)o < 1.25 mag. 

1.3. Aim of this paper 

The aim of this paper is to extend the /-band SBF calibration 
by Tonry et al. to blue colours and at the same time to test com- 
peting SBF models in that colour range. To this end we present 
SBF data for 25 Fornax cluster dwarf elliptical galaxies (dEs) 
in the range -16.5 < My < -11.2 mag, 0.8 < (V - /)o < 1.10 
mag. 

The paper is structured as follows: in Sect.|2lwe compare sev- 
eral SBF model predictions and correspondingly define four 
calibration test cases. Sect. |3] presents the technical details of 
the SBF measurements. In Sect. 0] the four calibration test 
cases are compared to the real SBF data and a purely empir- 
ical calibration is derived. Sect. [S] contains the discussion. We 
finish this paper with the conclusions in Sect.|6l 



2. SBF models 

SBF measurements of low surface brightness dEs have rather 
high error ranges, both in colour and in SBF magnitude, see 
Sect. l3.2l and Mieske et al. (|2003ai. Because of this we do not 
only aim at an empirical SBF calibration, but will also test ex- 
isting SBF model predictions with our measured SBF ampli- 
tudes, i.e. a semi-empirical approach. 

In Fig.[n theoretical values for M/ are plotted vs. (V - I)o for 
a set of old to intermediate age stellar populations with a wide 
range of metallicities, taken from diff'erent literature sources: 
Worthey (1994i, who uses his own stellar population mod- 
els based on VandenBerg (.1985 1 and Yale (Green et al. ll987> 
isochrones; Liu et al. J2000> . who use the updated models by 
Bruzual & Chariot (I993 l plus the BertelU et al. (T994\ Padova 
isochrones; Blakeslee et al. (2001j who use the stellar popula- 
tion models of Vazdekis et al. (I1996> and the updated Girardi 
et al. (2000 1 Padova isochrones; and Cantiello et al. ( I2003t . 
who use the stellar population synthesis code by Brocato et 
al. J 19991 12000> and the Teramo-Pisa-Rome isochrones (e.g. 
Castellani et al. 119911 Castellani et al. 119921 and Bono et 
al. I1997al I1997b> . Equation (□ is also shown in Fig. □ The 
Worthey model predictions are spHt up into the original predic- 
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tions from Worthey ( II 9941 1 ' and those that use the alternate stel- 
lar evolutionary isochrone library from Bertelli et al. (1]_994)^. 
For the sake of clarity, we separate the model predictions into 
two different plots. 

For red colours, all models trace equation Q reasonably well, 
with the original Worthey-models and those of Blakeslee et 
al. being more deviant from the empirical calibration than the 
other models. In the blue range, there is a notable discrepancy 
between the model predictions: the original Worthey ( 1994j 
models and those of Liu et al. (2000) predict a colour inde- 
pendent Ml ^ -2.0 mag for (V - /)o < 1.10 mag. In con- 
trast, Worthey models using the Bertelli et al. (!1994'i Padova 
isochrones, and also the models of Blakeslee et al. (2001 1 and 
Cantiello et al. (|2p03 1 show a large spread in M/ caused by 
age differences. Those models predict a continuation of equa- 
tion (0 for intermediate ages and a flattening of the relation for 
very old ages. At (V - /)o - 0.90 this leads to an uncertainty of 
up to 0.5 mag in relating M/ to (V - /)o. 
Based on the different theoretical model predictions in Fig.^ 
we define four "test cases" for possible calibrations: 
Case A: This is the "original" Tonry calibration, equation (Q, 
applied to all galaxies. 

Case B: This is equation (Q applied for (V - /)o > 1.09, and 
the following equation (|2j - which has half the slope of equa- 
tion (HJ - appUed for (V - /)o < 1.09: 

Mi = -2.0 -h 2.25 x ((V - /)(, - 1 .09) mag (2) 

This case represents the theoretical SBF models plotted in the 
left panel of Fig. ^ which viewed as a whole suggest a non- 
zero but shallower slope than that of the empirical calibration 
at redder colours. See case C) below for the reason of choosing 
{V-I)-l .09 as colour limit. 

Case C: This is equation ([Q apphed for (V - /)o > 1.09 and 
a constant of M/ = -2.00 mag apphed for (V - /)o < 1.09. 
This represents the theoretical SBF models plotted in the right 
panel of Fig. ^ The specific limit of (V - /)o = 1.09 is chosen 
because at this colour, the colour independent M/ ^ -2.0 mag 
of the models in the right panel of Fig.^matches the value of 
Mi from equation (Q). This is also the reason for adopting the 
same limit in case B), since the latter case represents the "hy- 
brid" solution between case A) and case C). 
Case D: A flat and a steep calibration relation are valid in the 
same colour regime. Depending on the strength of the SBF sig- 
nal, one of these is applied to a given galaxy. This case rep- 
resents the "two branch" hypothesis put forward in the inves- 
tigations by Jerjen et al. The steep relation is equation (0. It 
is apphed to galaxies which have (V - /)o > 1.09 or whose 
distance in case B) is smaller than the mean distance. The flat 
relation is simply a constant M/ - -2.00 mag, like in case 
C). This is applied to all galaxies with (V - /)o < 1.09 whose 
distance in case B) is larger than the mean distance. In other 
words, galaxies with m/ above the line for case B) in Fig. [2 
are assigned M/ = -2.00 mag, while those below are assigned 
relation (^. The mean distance obtained with this case is of 
course very close to that in case B). 

' http://astro.wsu.edu/worthey/dial/diaLa_model.html 
^ http://astro.wsu.edu/worthey/dial/diaLa_pad.html 



Those four cases will be tested against the real SBF measure- 
ments in Sect.|3] 

3. The SBF data 

The observations for this paper were taken in the nights 25th 
of October 2003 and 5-7 of December 2004 at Las Campanas 
Observatory (LCO), Chile. The instrument used was the 
Inamori Magellan Areal Camera and Spectrograph "IMACS" 
in imaging mode with the "short" f/2 camera, mounted at the 
6.5m Baade telescope. The double-asphere, glass-and-oil-lens 
f/2 camera produces an image of 27.4' field diameter at 0.20 
arcsec per pixel. The field is vignetted in the corners (by the 
tertiary mirror and its mounting assembly), changing from 0% 
flux loss at R = 12' to 10% at R = 15'. The IMACS detector 
is an 8192x8192 CCD mosaic camera which uses 8 thinned 
2kx4k detectors. Gaps of about 50 pixels separate the chips. 
In total, 7 fields in the central Fornax cluster were observed in 
the two bands V and /. The total integration times in V were 
between 1200 and 1800 seconds except for the central field 
around NGC 1399, for which it was 900 seconds. The total 
/-band integration time was between 3800 and 5160 seconds, 
except for the central NGC 1399 field, which had 1800 seconds 
integration. 

Standard star images of the Landolt fields SA92, SA95, SA98 
and SA104 - the latter one only in the second run - were taken 
at different airmass values in one short (2 seconds) and one long 
exposure (10 seconds). This resulted in about 110 data points 
for the first run and about 140 for the second one, each time 
corresponding to about 40 standard stars. Comparison of cali- 
brated colours and magnitudes for objects imaged in both runs 
yielded absolute zero point uncertainties of the order 0.02 mag. 

3.1. Image reduction before SBF measurement 

The image reduction steps before the SBF measurement were 
the following: first, a master-bias was created for each chip 
and was subtracted from the domeflat exposures. Then for each 
chip the bias corrected dome flats were combined to master- 
domeflats. Having the master-biases and master-domeflats pre- 
pared for each chip, we used the COSMOS package^ to do 
bias-subtraction, trimming and flat-field correction of the raw 
science frames. The reduced single science frames were reg- 
istered with integer pixel shifts to avoid distortion of the SBF 
power spectrum and combined using a min-max rejection al- 
gorithm. The seeing FWHM of the combined images generally 
ranged between 0.6 and 0.9", except at the edges of the field of 
view: especially for the first run in 2003, when the atmospheric 
dispersion corrector had not been installed yet, distortions close 
to the image border caused seeing degradations to about 2". 
There were two Fornax dEs for which we could not measure 
SBF because they were located in these regions of high distor- 
tion. 

The same image reduction procedure as for the single science 
frames was performed for the standard star images. Their in- 

^ Carnegie Observatories System for MultiObject Spectroscopy, 
|Ettp://llama.lco.cl/aoemler/COSMOS.html_ 
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Fig. 1. Theoretical "isochrones" in the M/-(y - /)o diagram from 5 different sources. For the sake of clarity, those were sep- 
arated into two different plots. Left: Blue short dashed lines from Worthey (1994 1 using Padova isochrones from Bertelli et 
al. .1994 (http://astro.wsu.edu/worthev/dial/dial_a_pad.htmli. Plotted are (from top to bottom in the blue range): 17, 12 and 8 



Gyrs. Metallicities range from -1.7 to dex from blue to red colours. Green long dashed lines from Blakeslee et al. (1200 1> 
for ages of 7.9, 12.6 and 17.8 Gyrs and metallicities from -1.7 to -1-0.2 dex. Magenta solid lines from Cantiello et al. (20031 
for ages of 5, 11 and 15 Gyrs and metallicities from -2.3 dex to 0.3 dex. Right: Blue short dashed lines from Worthey U994t 
http://astro.wsu.edu/worthey/dial/dial_a_model.htmli. Plotted are 17, 12 and 8 Gyrs. Metallicities range from -1.7 to dex from 
blue to red colours. Red long dashed lines from Liu et al. (2000 1 for ages 8, 12 and 17 Gyrs and metallicities between -2.3 and 
0.4 dex. In both panels, the solid line for {V - I)o > 1 .09 indicates the empirical calibration (Q by Tonry et al. ( I1997> . the solid 
line for {V - I)o < 1.09 indicates equation (|2j, whose slope is half that of equation (Q. The dashed lines indicate a constant 
Ml = -2.00 mag and the continuation of equation ([0 for (V - /)o < 1.09. See Sect.l^for more details. 



strumental magnitudes were measured with the IRAF package 
APPHOT in apertures equal to those used by Landolt (1994i. 
Then, a single photometric solution was determined for the 
entire 8kx8k image. Chip-to-chip variations around the mean 
photometric zero point were smaller than 0.02 mag. To cor- 
rect for galactic reddening and absorption, we used the val- 
ues from Schlegel et al. (1998 1, who give A/ = 0.025 and 
E{V - /) = 0.018 for the coordinates of the Fornax cluster 

3.2. SBF measurement 

The SBF measurement procedure was similar to that already 
outlined in Mieske et al. (.2005J based on VLT FORS 1 photom- 
etry. Some nuances have been modified: the fluctuations from 
the sky background were now derived in the same image close 
to the investigated dE instead of in a comparison field. Those 
sky fluctuations consisted of undetected background galaxies 
and fringing. The seeing was not constant over one chip, gen- 
erally increasing towards the edges of the field of view. This 
caused a correspondingly varying amount of fluctuation from 
undetected background galaxies along the chip. Also the fring- 
ing varied in some cases significantly over one chip, see Fig.|3 
These two dependences on position required a local determina- 



tion of sky fluctuations. Note that for the VLT FORS photom- 
etry there was no significant fringing nor strong variations of 
seeing among the different exposures. 

In turn we describe the relevant measurement steps for deter- 
mining the SBF magnitude: 

1 . Model mean galaxy light with ELLIPSE using a sigma clip- 
ping algorithm to disregard contaminating sources, subtract the 
model. 

2. Detect and subtract remaining contaminating objects from 
original image. 

3. Model mean galaxy light on the cleaned image. 

4. Subtract model of original image. 

5. Divide resulting image by square root of the model, cut out 
circular portion with radius typically 20-25 pixel (4 - 5"), cor- 
responding to about 4-8 seeing disk diameters. 

6. Mask out contaminating sources like foreground stars and 
background galaxies. 

7. Calculate the power spectrum (PS) of the cleaned image. 

8. Calculate the PS of the sky background in three different 
blank image sections close to the investigated galaxy. This 
sky background PS contains both fluctuation contributions 
from undetected background galaxies as fluctuations caused 
by fringing, see explanation at beginning of this subsection. 
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Normalise these fluctuations by dividing these blank images 
by the mean galaxy intensity in the region where SBF are mea- 
sured. Subtract this normalised background PS BG from the PS 
of the SBF image. 

9. Obtain the azimuthal average of the resulting PS. 

10. Fit function of the form 

P(k) = E{k) xPo+Pi (3) 

to die result. Here, E{k) = PSF(k) » W{k) with W{k) the PS 
of the mask used to excise contaminating sources and restrict 
the measurement area. PS F{k) is the undistorted power spec- 
trum of the PSF normalised to unity at k=0. PSF{k) is derived 
from one to three well exposed star(s) in the same image close 
to the respective galaxy. When more than one star was used, 
the average PSF-PS was adopted as PSF(k). The variations in 
PSF FWHM between stars chosen for one galaxy were small, 
around 5%. Using E(k) instead of PS F{k) implicitly corrects 
for the power spectrum damping at low wavenumbers that is 
caused by the convolution of the undistorted PS F(k) with the 
power spectrum of the mask (see Mieske et al. ^005). The 
amount of power spectrum damping at wavenumber zero was 
on average 0.11 mag with a galaxy-to-galaxy dispersion of 0.04 
mag. The signal-to-noise ratio of the SBF measurement is then 
defined as S /N - ^. Since the background fluctuations have 
been subtracted from the power spectrum before the fitting, this 
definition corresponds to the ratio of pure SBF signal over sky 
noise. For the SBF analysis, only galaxies with S /N > 3 were 
considered (see for example Figs|3]and|5j- 

1 1 . Obtain the desired observable m/ from 

m, = -2.5 * logiPo) + ZP, -A,-Ak + Aqc (4) 

A/ is the foreground absorption, = z x 7 is the k-correction 
for SBF in the /-band (Tonry et al. 19971. For calculating Ak 
we assumed cz - 1500 km/s for all galaxies. Agc is the fluc- 
tuation contribution from undetected globular clusters (GCs). 
To estimate Agc, we determined the completeness magnitude 
for detecting unresolved point sources for each galaxy sepa- 
rately according to the seeing via artificial star experiments. 
The 50% completeness magnitudes ranged between 1-22.5 and 
24.2 mag. For calculating Agc, we assumed equation (15) from 
Blakeslee & Tonry (I1995> . a turn-over magnitude (TOM) of 
Mj - -8.46 mag and Gaussian width cr = 1.2 mag for the 
GCLF (Kundu & Whitmore l2001> . a Fornax distance modu- 
lus of 31.4 mag, and a specific frequency of its GC system 
Sn ,o-„S,,i5) = 5 (MiUer et al. 1998). The calculated val- 
ues for Agc are listed in Tabled ranging between 0.01 mag and 
0.42 mag with a mean of 0.06 mag and median of 0.03 mag. 
The colour (V - /)o of each galaxy was adopted as the mean of 
two values: first, the {V-I)o estimates from the IMACS data in 
the region where SBF were measured; second, the colours de- 
rived by us for the same galaxies in Hilker et al. ( I2003I I. These 
latter data were based on WFCCD photometry with the 2.5m du 
Pont telescope, also at Las Campanas. Those colours had been 
determined in circular apertures of 4" radius, which is in close 
agreement to the size of the region where SBF were measured. 
For those galaxies entering our SBF data set, the systematic dif- 
ference between the IMACS and WFCCD colours is 0.017 + 



0.021 mag, with a scatter of 0.10 mag. This shows that our pho- 
tometric colour calibration is accurate to better than 0.02 mag. 
Fig-Eshows the CMD of the investigated galaxies for the var- 
ious colour sets. The rms dispersion in colour of a linear fit is 
0.083 mag for the IMACS data, 0.068 for the WFCCD data and 
0.058 mag for the mean, independent of magnitude. We adopt 
this very dispersion of 0.058 mag as colour uncertainty. This is 
to be seen as an upper limit to the real colour error because the 
rms dispersion certainly includes some "cosmic" scatter 
With (V - /)o at hand, Mj was derived using equation Q, 
which is calibration case A). Together with the measured m/, 
this yielded the "null hypothesis" distance modulus (m - M)a, 
see Table[n Fig.0]shows thumbnail images illustrating the SBF 
measurement procedure for two out of the 28 dEs with clear 
SBF signal. 25 of the 28 investigated galaxies hadS /N > 3 and 
entered our subsequent analysis. Their absolute brightnesses 
are in the range - 16.5 < My < - 1 1 .2 mag, their central surface 
brightnesses cover the regime 20.3 < /uy < 24.6 mag/arcsec^. 
They occupy a colour range 0.8 < (V - I)q < 1.10 mag, ide- 
ally suited for extending the Tonry colour range to the blue. In 
Fig.|5]we show that the null hypothesis SBF distance modulus 
does not correlate with the S/N of the SBF measurement nor 
with the amount of background fluctuations BG. 
Note the large range of BG values in Fig.|5] This has three rea- 
sons: The first and main reason is related to the fact that BG is 
determined by fitting equation ^ to the normalised sky back- 
ground PS. For this fitting we have to exclude the same low 
wavenumber regime that is excluded for fitting the galaxy SBF 
amplitude Pq. Since BG is only a fraction of Pq, this exclusion 
of low wavenumbers has a stronger effect on the fitting preci- 
sion for BG than for fitting Pq- Indeed, the negative value for 
BG that is measured for FCC 215 (see Tableland Fig.O is 
a consequence of this low wavenumber exclusion. The back- 
ground PS in the remaining wavenumbers is predominantly 
noise dominated such that a formally negative value for BG is 
fitted. As a consequence of that, the SBF data with poor seeing 
(FWHM>0.85") show a scatter in BG of 0. 16 mag which is al- 
most twice as large as the scatter of 0.09 mag for the data with 
good seeing (FWHM<0.85"), see Fig.|5] The second reason 
for the variation of BG is that the galaxies investigated span a 
considerable range in surface brightness (> 4 mag). This leads 
to varying relative amounts of BG, given that the SBF signal in 
lower surface brightness regions is more strongly affected by 
sky background fluctuations. The third reason is the somewhat 
varying fringing amplitude, both between different chips and 
different nights, see for example Fig.|2 

We finally note that the sky background correction of the 
galaxy SBF data is done by directly subtracting the sky back- 
ground PS image from the normalised galaxy light PS image. 
Only the quantification of the correction is achieved hy fitting a 
PS. To take the uncertainty arising from the power-spectrum fit 
and from the sky fluctuation variability into account, we adopt 
the scatter of Pq resulting from the subtraction of three differ- 
ent sky background PS images as error estimate of the SBF 
measurement, see the following Sect. 13. 2.11 
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Fig. 2. Two image excerpts from the same chip and exposure close to two investigated dEs, illustrating the variable degree of 
fringing. In the left thumbnail, FCC 194 is at the bottom right. In the right thumbnail, FCC 196 is at the bottom. Both images 
have the same intensity cuts. Note the fringing patterns in the right thumbnail, which are not present in the left one. 
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Fig. 3. Colour magnitude diagram of the galaxies entering our 
SBF analysis, i.e. those ones with an SBF S/N > 3, see also 
Fig. |5l Small blue circles indicate colours derived from the 
present IMACS data set. Small red asterisks indicate colours 
from the WFCCD data set of Hilker et al. |2003 l. Large encir- 
cled asterisks is the mean of both values. The V magnitude is 
from the WFCCD data set for all galaxies. 



The error of Pq was estimated in two ways: 

1 . based on the Monte Carlo simulations presented in Mieske 
et al. (2003a I. For the median seeing of our data (0.75"), inter- 
polation of the simulation results for 0.5 and 1.0" for Fornax 
dEs and 1 hour integration yields the following uncertainties: 
0.37 mag in the range -11 > My > -13 mag, 0.27 mag in 
the range -13 > My > -14.5 mag, and 0.20 mag in the range 
My < -14.5 mag. Note that the simulation errors are calcu- 
lated for images of the same pixel scale than the one used here. 
The fainter VLT-FORS zero points assumed in the simulations 
are practically compensated by the larger integration times of 
our Magellan-IMACS data. 

2. from the scatter of the fitted values for Pq when subtracting 
the three different sky background power spectra. This was an 
important double-check of the errors adopted from the simu- 
lations, because the fringing of the IMACS data represent an 
additional uncertainty source, apart from the photon noise and 
background galaxy fluctuations which had been included in the 
simulations. 

The maximum of both error estimates was adopted as error in 
Po- For seven out of the 25 dEs investigated, the scatter from 
background subtraction was larger than the Monte Carlo esti- 
mate, see Table[2 The resulting mean error of (m - M) is 0.41 
mag when adopting calibration case A). There are four galax- 
ies in the range -16.5 < My < -11.5 mag which were im- 
aged twice in the course of our survey since they were located 
in the overlap of adjacent fields. The differences in (m - M) 
between those double measurements ranged between 0.10 and 
0.35 mag. This is encouragingly low, indicating that our uncer- 
tainty estimates are on the conservative side. 



3.2.1 . Error estimate for SBF distances 

The error of (m - M) has two equally important contributions: 
uncertainty in (V - /)o and Pq. 

The uncertainty in (V - 1\) of 0.058 mag (see above) translates 
into a 0.26 mag distance uncertainty in the case of the empiri- 
caUy calibrated slope 4.5. 



3.3. Correlation between seeing and mi 

A final test for observational biases before coming to the cal- 
ibration is to look at the measured SBF amplitude Irij as a 
function of seeing, see Fig. |6l Those two observables should 
be independent of each other However, we do find a 2cr sig- 
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Table 1. Results of SBF measurements for the 25 Fornax dEs with S/N> 3 in the SBF measurement. The galaxies are ordered by 
decreasing brightness. In the first column, the reference number in the Fornax Cluster Catalog FCC (Ferguson & Sandage [l989> 
is given. The last four galaxies had been too faint for this catalog. They are labelled with a field number and running index 
according to Hilker et al. (2003 1. All columns have magnitudes as units except S/N and seeing. (V - /)o is the mean colour of 
the photometry from Hilker et al. J2003> and this paper. V is from Hilker et al. J2003t . The other parameters are derived in this 
paper. BG gives the amount of background fluctuation present in the original fluctuation image, comprising both fluctuations 
from undetected background sources and fluctuations arising from fringing, (m - M)^ to (m - M)^ give the four sets of distances 
obtained when applying the four different calibration test cases, see Sect.|2 



nificant coiTelation in the sense that bad seeing data have 
stronger m/. We neither find a coiTelation between seeing and 
the background fluctuation amplitude BG nor between seeing 
and galaxy luminosity. 

One possibility to qualitatively explain the trend of seeing with 
nil is that a fit to the same seeing power spectrum yields sys- 
tematically higher Pq when excluding more low k data points 
(which is the case for bad seeing data). We have tested this 
for our stellar power spectra: for a fixed assumed width of the 
power spectrum, the fitted Po depends only very marginally (to 
a few percent) on the limiting fit wavenumber k, even when the 
fitted power spectrum data points are all below 1/3 of Pq- We 
have also double-checked whether for some large seeing data 
we may erroneously have used marginally resolved objects for 
PSF fitting. However, this was not the case. 
Jensen, Tonry & Luppino (1998 1 argue that poor quality SBF 
data in IR SBF measurements are prone to artificial increases of 
the fluctuation amplitude due to sky background fluctuations, 
imperfect galaxy modelling or undetected background sources. 
For our data, we can most likely exclude the effect of sky fluc- 
tuations and undetected background sources, see Sect. 13. 21 The 
sky fluctuations are independent of seeing and derived distance 
(see above), and the globular cluster fluctuations are small. We 



conclude that the imperfect galaxy subtraction (see e.g. right 
panel of Fig.|^ is the major contributor to the increase of fluc- 
tuation power at large seeing. This is because the PSF scale 
approaches that expected for galaxy model residuals (the SBF 
images for the worst seeing data have a diameter of only about 
8 PSF FHWM) such that the commonly applied rejection of 
low wavenumbers in the power spectrum fit can apparently not 
remove all the residual power for the bad seeing data. 
For the SBF calibration, we therefore exclude all data with see- 
ing worse than 0.85" FWHM. Doing this removes the depen- 
dence of seeing on m/, see Fig. |6l This exclusion applies to 
eight galaxies, such that we are left with a final sample of 17 
SBF data points for the calibration. In the remainder of the pa- 
per the seeing restricted calibration is the reference one, but we 
will also indicate the calibration results for the full sample of 
galaxies. 

4. SBF calibration 

AA. The semi-empirical approach 

For each of the four calibration test cases defined in Sect.|2 the 
corresponding distribution of SBF distance modulus vs. V and 
(y - /)o is shown in Figs.0and|8l Fig.|9lplots m/ vs. (V - /)o. 
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Table 2. Mean SBF distance moduli for the 17 investigated dEs with seeing better than 0.85", assuming four different calibration 
cases (Sect. |2j. cr^m-M) gives the scatter of the measured distances around their mean. 6(„,-m) gives the mean uncertainty of a 
single distance measurement. The two last columns show the correlation between (m - M) and (V - I)o, and (;« - M) and Vo, 
see also Figs.0and|8l The last line gives the Fornax cluster reference distance from the HST-Key project (Freedman et al. l2001> . 
Note that all errors indicated are 1 cr except the eiTor range for cr^m-M)^ which is a 90% confidence range. 
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Table 3. Like Table|2] but now including all 25 investigated dEs. 



In all plots, galaxies with seeing worse than 0.85" are indi- 
cated. Table[ngives the results in numbers. Table|2lsummarises 
the parameters for each calibration case when using only those 
17 galaxies with seeing better than 0.85". Shown is the mean 
(m - M), the distance scatter, the mean measurement error, and 
the dependence of {m - M) on colour and magnitude ( 

and ^^f^). TableEJshows the same, but for all 25 galaxies. 
As a starting point to test these cases, we have to adopt a fidu- 
cial Fornax cluster distance modulus. The mean of the Cepheid, 
PNLF, GCLF and /-band SBF distance from the HST-Key 
project summary of Ferrarese et al. (2000( is (m - M) - 31.52 
mag. Note that here, the Cepheid scale is used to calibrate 
the other methods. Freedman et al. (2001 1 present a revised 
Cepheid calibration, which causes a downward correction of 
0.13 mag for the average Cepheid distance to the Fornax, 
Virgo and Leo 1 clusters (Table 8 of their paper). The result- 
ing corrected mean Fornax distance then is 31.39 mag. This 
corresponds exactly to the revised Cepheid distance for Fornax 
from Freedman et al. (2001 1, which has an error of 0.12 mag. 
Therefore, we define a reference distance of {m - M) - 31.39 + 
0.12 mag. Because of the same Cepheid scale revision, also the 
formal zero-point of equation Q becomes fainter by 0.09 + 
0.02 mag, as calculated from comparing the revised Cepheid 
distances with the old ones used by Tonry et al. ( 1997 1 for their 
group calibration. However, this is a small change compared to 
the zero-point scatter of models shown in Fig.^ We did there- 
fore not include this offset in the test cases, also because the 
Tonry calibration did not extend to blue colours. 
We define three control criteria for the four different sets of dis- 
tance moduli: 1) we demand that there is no significant correla- 
tion between colour (V-I\) and distance {m-M). Since we aim 
at a calibration of Mj with (V - I)o, this is a very fundamental 
requirement. 2) we require that the scatter of the distance val- 



ues around their mean is not significantly smaller than the mean 
single measurement uncertainty. This is to avoid that we over- 
correct possible trends in our data. 3) we demand that there is 
no significant correlation between Vq and (m - M). Due to the 
well known colour-magnitude correlation for galaxies (see also 
Fig.|3}, this is an important double-check of criterion 1). 
Of the four test cases, the two cases B) and D) agree best with 
the reference distance. The null hypothesis case A) yields a 
marginally higher distance (1.4cr), the "constant" case C) a 
marginally lower one (1.7cr). The difference between measure- 
ment scatter and mean single measurement error is largest for 
case C). This is a 2.2cr effect, consistent with non-negligible in- 
trinsic scatter of the caUbration ("cosmic variance"). As men- 
tioned before, we are confident that our estimate of the sin- 
gle measurement uncertainty is not too optimistic. Therefore 
we take the disagreement for case C) as a hint that a con- 
stant Mi with small cosmic variance - as predicted by the Liu 
and Worthey models, see Fig.^- does not fit our SBF mea- 
surements. The large scatter of case C) is accompanied by a 
clear trend of V with (m - M), significant at the 3cr level. This 
strengthens the conclusion that case C) is not the best approx- 
imation of our data. Also in case B) there is some correlation 
between V and (m - M), but at slightly lower significance (2cr). 
Case A) shows a marginal trend of colour with distance, but 
only at the 1.3cr significance. For the shallower calibration 
cases B) and D) the nominal value of the trend becomes very 
small. For case C) the nominal slope is comparable but of oppo- 
site sign to that in case A). This is consistent with a slight flat- 
tening of the relation between M; and (V - /)o at bluer colours. 
Of the four test cases, case D) shows the best agreement with 
the data: the highest disagreement is a very marginal 1 . Icr trend 
of (m - M) with V. The low distance scatter in case D) is a nat- 
ural consequence of the adoption of two separate calibration 
branches over a common colour range. 
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Fig. 4. Example thumbnails of the SBF measurement proce- 
dure for two Fornax cluster dEs. Left FCC 194, right LSB6-4. 
Thumbnail size is 12.5". Image sequence from top to bottom: 
1. Original galaxy image. 2. Galaxy light model subtracted 
from former image. 3. Former image divided by square root of 
galaxy model, with contaminating sources/fit residuals masked 
and with the measurement region restricted to a circular area. 
4. Azimuthally averaged power spectrum of the former SBF 
image, with the fit indicated by the dashed line. 
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Fig. 5. Distance modulus for case A) (see text) plotted vs. 
S/N of the SBF measurement (left) and the amount of back- 
ground fluctuations BG contributing to the measured SBF sig- 
nal (right). The three sources marked with large circles in the 
right panel are those below S IN - 3 in the left panel. Filled 
circles are for galaxies with seeing below 0.85", open circles 
for seeing above 0.85", see Fig.|6land text. 
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Fig. 6. The fluctuation magnitude m/ plotted vs. the seeing 
FWHM of stars in the vicinity of the dEs investigated. There 
is a 2cr correlation between both entities, in the sense that bad 
seeing images are biased towards too bright SBF magnitudes. 
For the SBF calibration, we use only measurement with seeing 
better than 0.85", as indicated by the dashed line. For those 17 
galaxies, the correlation between seeing and m/ is insignificant. 



We note that for the full sample, case C) is rejected at even 
higher significance due to a lower mean distance and stronger 
correlation between (m - M) and (V - /). 



4.2. The empirical approach 

We also give the best fit result of a purely empirical calibration 
valid in the colour range 0.85 < (y-/)o < 1.10 mag, assuming 
(m-M) = 31.39 + 0.12 mag: 

M/ = -2.13(+0.17) + 2.44(+1.94) x [(V - /)o - 1.00] mag (5) 
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Fig. 7. Apparent V magnitude plotted vs. SBF distance modu- 
lus for the four different assumed calibration relations A) to D), 
as defined in Sect.|2]and Table|2] Filled circles are galaxies with 
seeing better than 0.85". Their mean distance and its Icr eiTor 
are indicated by the vertical solid and dashed lines. The green 
ticks at (m - M) - 31.39 mag mark the reference Fornax cluster 
distance from Freedman et al. ( 2001 ) and its error. Open circles 
are data points for galaxies with seeing worse than 0.85". They 
were not included in calculating the mean distance and are only 
shown for the sake of completeness. 
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Fig. 9. Apparent SBF magnitude m/ plotted vs. (V-l)o colour 
for the same galaxies as in Figs. and |S1 Overplotted are the 
calibration relations from Fig.[2for a distance modulus of (m - 
M) - 31.38 mag, which is the one obtained in case B). Solid 
line plus dotted line is calibration case B). Case D) is defined 
by applying the solid and short dashed line to galaxies with 
(y - /)() < 1.09 and ot/ fainter than the line for case B), while 
applying the sohd and long dashed line to the other galaxies. 
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Fig. 8. (V-I)o colour plotted vs. SBF distance modulus for four 
different assumed cahbration relations, see Fig.Q 



The zero point error includes our statistical errors and the error 
of the Fornax cluster reference distance. The statistical errors 
of the zero point and the slope were derived from the disper- 
sion of the data points around the linear relation, not from the 
single measurement eiTors (this would have decreased the er- 
rors by 50%). The formal disagreement with the Tonry slope of 
4.5 is l.lcr, while the disagreement with a slope is 1.3cr. For 
the full sample of galaxies, the zero-point changes to -2.31 + 
0.15, and the slope to 2.80 ± 1.36. Not surprisingly, this agrees 
with the finding from the semi-empirical approach that case C) 
is rejected at higher significance for the full sample. 
To quantify the intrinsic scatter - the "cosmic variance" - of the 
empirical "good seeing" calibration, we subtract in quadrature 
the mean single measurement uncertainty in m/ of 0.30 mag 
from the dispersion around the fit, which is 0.46 mag. The re- 
sulting cosmic variance is 0.34 + 0.14 mag. This is larger than 
the value of 0.05-0.10 mag estimated by Tonry et al. J1997> 
in their survey at redder colours, in agreement to those models 
that predict M/ to be more sensitive to different age/metallicity 
combinations in the blue than in the red (see left panel of 
Fig. [Q. Note that the two theoretical predictions by Liu et al. 
and Worthey et al. in the right panel of Fig. ^that predict a 
colour independent M/ (case C)), do not exhibit a significant 
scatter. 
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5. Discussion 

The results of the previous section suggest an overall shallower 
slope than the Tonry slope determined for red colours, accom- 
panied by a significant cosmic scatter. The question is whether 
this is created by a single relation with some scatter as in case 
B) (see also Mei et al. 2005 1, or by a two-branch solution as in 
case D) (Jerjen et al. l2001ll l. Case D) contains the implicit as- 
sumption that a) the relation between M/ and (V -I)o broadens 
up significantly at blue colours; b) the age-metallicity distribu- 
tion of the galaxies investigated is broad; and c) the investigated 
galaxies consist of two distinct populations: an old, metal-poor 
one plus an intermediate age, more metal-rich one. Such a bi- 
modality in stellar populations is known from the colour dis- 
tribution of globular cluster systems (e.g. Peng et al. 2005' 
Kissler-Patig 2000 Harris et al. 2005 Dirsch et al. 2003 1, but it 
has not been found for bright dEs (e.g. Geha et al. 2003 Held 
& Mould 1994 1. In order to support a two-branch calibration in 
the sense of the Jerjen papers - as opposed to an overall shal- 
lower relation with some intrinsic scatter - there have to be two 
well separated galaxy populations in the Mj-(V - /)o plane. 
To investigate this, we test the distance distribution in the dif- 
ferent calibration cases for bimodality. We use KMM (e.g. 
Ashman, Bird, & Zepf 1994) in homoscedastic mode to quan- 
tify the probability with which a double-peaked Gaussian dis- 
tance distribution is preferred over a single-peaked one. This 
is certainly a very crude estimator, but the measurement uncer- 
tainties and low numbers do not allow for more detailed consid- 
erations. The confidence levels with which a bimodal Gaussian 
is preferred over a unimodal one are the following: 88% for 
case A), 92% for case B), 64% for case C), and 27% for case 
D). The strongest case for a double-peaked distribution is case 
B), yielding one main peak with 12 galaxies and a secondary 
one with five. These five galaxies are the ones below the long 
dashed line in Fig.|9l The fact that the probabilities for cases B) 
and D) are so different is consistent with a proper two-branch 
calibration as opposed to a unimodal one with broad scatter. 
However, the confidence level for case B) is only 1.8cr, and 
therefore this finding is still of indicative nature. We note that 
the SBF S/N and also the amount of background fluctuations 
BG is not a function of distance modulus. Both entities are in- 
distinguishable between SBF-bright and SBF-faint sample. 
Are there any other diff'erences between dEs belonging to the 
bright and faint m/ part of case D)? The difference in mean 
colour (V - I) between SBF-faint and bright dEs is insignifi- 
cant, being -0.015 + 0.028 mag. However, the mean V-band to- 
tal luminosity of the SBF-faint sample is 1 .6 ± 0.5 mag brighter 
than that of the SBF-bright sample, accompanied by a 1 mag 
brighter mean central surface brightness and 40% larger expo- 
nential scale length. An obvious conclusion from this finding 
is that the lower luminosity and smaller dEs in Fornax have in 
comparison to their brighter counterparts a younger integrated 
age and somewhat higher metallicities, see Fig. ^ This may 
imply that the Fornax dwarf galaxy populations were built up 
in subsequent events, where the most recent events produced 
less massive dwarf galaxies than at earlier times. The 92% sig- 
nificant bimodality in m/ is in that context consistent with at 
least two separate dwarf formation episodes in Fornax. 



We finally note that the average projected radial distance to 
NGC 1399 is indistinguishable for both samples, indicating 
that both groups are not distributed in a significantly differ- 
ent fashion. Due to the intrinsic faintness of the investigated 
dEs, only three of them have spectroscopic metallicities avail- 
able (Held & Mould [T994I and Mieske et al.|2^6 i, of which 
only FCC 211 is bluer flian (V - I)o = 1.09. This obviously 
does not allow one a spectroscopic confirmation of bimodality 
in the age-metallicity distribution. 

5.1. Revisiting the Hydra and Centaurus SBF 
measurements 

At the end of this discussion we revisit the relative distance be- 
tween Hydra and Centaurus. The purely empirical calibration 
relation (|5) for 0.85 < (V - /)o < 1.10 mag can be used to 
check the SBF distances to the dE dominated sample of Hydra 
and Centaurus cluster galaxies from Mieske & HiUcer ('2003b^ 
and Mieske et al. (2005 1. There is a 0.05 mag overlap of the 
colour regime for relation ^ with that of Tonry's relation Q. 
In the Hydra/Centaurus papers we assumed relation Q to hold 
for (V-I)o > 1.00, and a 50% smaUer slope for (V-I)o < 1.00. 
Both the zero-point and slope of relation (|5} are consistent with 
that adopted for {V-I)o < 1 .00 in those papers. That is, the SBF 
distances remain unchanged in the case of {V-I)o - 1 .00 as the 
limit between steep and shallow slope. If we restrict the validity 
of relation Q to Hydra/Centaurus galaxies with (V-I)o > 1 . 10 
and adopt the slope of relation (|5} for bluer colours, the mean 
distance of both samples becomes 0.09 mag lower. I.e. there is 
no consequence for the relative distance. The uncertainty of re- 
lation's (|5j zero-point is larger than this difference. Therefore 
we can state that the blue SBF calibration from this paper does 
not change the relative Hydra/Centaurus distance nor does it 
require a change in the absolute distances. 

6. Summary and conclusions 

We have presented /-band SBF measurements taken with 
IM ACS ©Magellan for a sample of 25 dEs in the Fornax clus- 
ter in the colour magnitude range -16.5 < My < -11.2 mag 
and 0.8 < (V - /)o < 1.10 mag. The aim of this investigation 
was to address the SBF calibration at blue colours. Our colour 
calibration is accurate to better than 0.02 mag. To avoid obser- 
vational biases towards too bright SBF amplitudes, we restrict 
our calibration sample to those 17 galaxies with seeing better 
than 0.85". 
These are our results: 

- Our SBF data are inconsistent at the 3cr confidence level 
with a constant, colour independent absolute SBF magni- 
tude Mi at blue colours, as suggested for example in the 
SBF models by Liu et al. (2000) and Worthey et al.jg593|l. 

- Our SBF data support a non-zero slope between M/ and 
(y-/)o. There is some evidence that this slope is shallower 
than the value of 4.5 found by Tonry et al. J1997> at redder 
colours, but at less than 2cr. 

- In the Mi - (V - I) plane, our SBF data favour at the 1.8 
cr confidence level a two-branch calibration relation over 
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a single one with some uniform scatter We find evidence 
that the fainter galaxies in our samples exhibit younger inte- 
grated ages and higher metallicities than the brighter galax- 
ies. 

- From a purely empirical SBF calibration (relation|5} we de- 
duce a 0.34 + 0. 14 mag cosmic scatter of M/. This is several 
times larger than found by Tonry et al. (1997i for redder 
colours, in agreement with those theoretical predictions that 
predict M/ to be more sensitive to different age/metallicity 
combinations in blue colours than at red colours. 

- Applying the empirical calibration to the SBF measure- 
ments of Hydra and Centaurus cluster dEs from Mieske & 
HiUser (|2b03bi and Mieske et al. (2005) has no effect on 
the relative SBF distance of both galaxy clusters. Also the 
absolute distances do not require correction. 
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